State-defect constraint pairing graph coarsening method for Karush-Kuhn-Tucker matrices arising in orthogonal collocation methods for optimal control

A state-defect constraint pairing graph coarsening method is described for improving computational efficiency during the numerical factorization of large sparse Karush–Kuhn–Tucker matrices that arise from the discretization of optimal control problems via an Legendre–Gauss–Radau orthogonal collocation method. The method takes advantage of the particular sparse structure of the Karush–Kuhn–Tucker matrix that arises from the orthogonal collocation method. The state-defect constraint pairing graph coarsening method pairs each component of the state with its corresponding defect constraint and forces paired rows to be adjacent in the reordered Karush–Kuhn– Tucker matrix. Aggregate state-defect constraint pairing results are presented using a wide variety of benchmarkoptimal control problemswhere it is found that the proposed state-defect constraint pairing graph coarsening method significantly reduces both the number of delayed pivots and the number of floating point operations and increases the computational efficiency by performing more floating point operations per unit time. It is then shown that the state-defect constraint pairing graph coarsening method is less effective on Karush–Kuhn–Tucker matrices arising from Legendre–Gauss– Radau collocationwhen the optimal control problemcontains state and control equality path constraints because such matrices may have delayed pivots that correspond to B Anil V. Rao [email protected] Begüm Şenses Cannataro [email protected] Timothy A. Davis [email protected] 1 Department of Mechanical and Aerospace Engineering, University of Florida, Gainesville, FL 32611-6250, USA 2 Department of Computer Science and Engineering Information and Science Engineering, Texas A&M University, College Station, TX 77843, USA